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ABSTRACT 
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1.   Introduction 

In  computational  geometry,  optimization  problems  are  often  posed  in  a  continuous 
(as  opposed  to  discrete  or  combinatorial)  setting.  One  can  resort  to  numerical  methods 
to  give  approximate  solutions  to  any  degree  of  accuracy  or  solve  the  problem 
symbolically  and  reduce  the  problem  to  root-finding  or  to  a  decision  procedure  for 
Tarski's  language  for  elementary  geometry  and  algebra  [BKR84].  A  recent  paper  by 
Sharir  and  Schorr  [SS84]  shows  a  case  where  no  combinatorial  finiteness  criterion  (except 
indirectly,  by  a  reduction  to  Tarslci's  language)  is  known:  it  is  the  problem  of  finding 
shortest  paths  between  a  pair  of  points  among  polyhedral  bodies  in  space.  The  preferred 
method  for  these  problems,  however,  is  to  find  a  combinatorial  "finiteness  criterion"  for 
each  problem.  One  example  where  such  an  approach  works  (extremely  well,  in  practice) 
is  linear  programming.  In  this  paper,  we  address  another  such  problem,  the  potato- 
peeling  problem  described  in  the  abstract.  This  problem  was  first  posed  (in  a  more 
general  form  than  we  have  stated  it)  by  J.  Goodman  [G008I]  who  obtained  various 
mathematical  properties  of  the  problem.  The  only  partial  solution  in  [G008I]  gives  the 
criteria  for  a  finite  solution  if  the  polygon  has  n  ^  5  sides.  Quite  independently,  T. 
Woo  [Woo83]  studied  the  same  problem  dubbing  it  the  "convex  skull"  problem.'  Not 
only  will  we  show  that  the  problem  is  finite  in  general,  but  we  will  derive  a  polynomial 
time  bound. 

Our  computational  model,  as  is  typical  in  this  subject,  assumes  a  random-access 
computer  with  infinite  precision  real  arithmetic.  Observe  that  the  notion  of  a  "finiteness 
criterion"  is  a  relative  one,  in  view  of  the  infinite  precision.  Furthermore,  we  assume 
that  the  solution  of  simple  trigonometric  equations  takes  0(l)  steps  (this  amounts  to 
assuming  the  availability  of  trigonometric  functions  and  their  inverses).    For  example. 


'  Indeed,  we  are  indebted  to  T.   Woo  who  first  brought  this  problem  to  our  attention.    Later,  M. 
Sharir  pointed  out  the  work  of  J.  Goodman. 


given  angles  a,  [3  and  constant  c,  we  assume  that  we  can  in  (9(1)  time  find  the  angle  8 
satisfying 

sin  (9  +  a)   ^ 
sin(e-p)        ^' 

It  turns  out  that  our  method  applies  to  other  problems  that  have  been  studied.  The 
general  framework  for  these  problems  can  be  posed  as  follows.  Let  P,  Q  be  families  of 
polygons,  and  let  jjl  be  a  real  function  on  polygons  with  the  property  that  for  all  P,  Q  in 

P  CQ     ^     MP)  ^  M-(Q)- 

Note  that  this  property  holds  if  jjl  measures  the  area.  If  the  polygons  in  Q  are 
convex,  the  property  also  holds  when  jjl  measures  the  perimeter.  In  this  paper,  polygons 
are  assumed  simple  in  the  sense  that  they  are  self-avoiding,  and  polygonal  regions  (also 
called  "polygons"  when  the  context  is  clear)  are  simple  in  the  sense  that  the  boundary  of 
each  region  forms  a  simple  polygon.  The  class  of  (polygon)  inclusion  and  enclosure 
problems^  are  defined  as  follows: 

InciP,  Q.,  \i.)  :  Given  P  i  P,  find  the  ^.-largest  Q  ^  Q.  that  is  included  in  P. 

Enc{P,  Q.,  \l)  :  Given  P  ^  P,  find  the  ^.-smallest  Q  ^  Q.  that  encloses  P. 

These  two  classes  of  problems  are  "duals"  in  some  sense  but  we  know  of  no 
systematic  way  whereby  an  algorithm  for  a  problem  can  be  transformed  to  one  for  its 
dual.  For  instance,  the  dual  of  potato-peeling  problem  is  the  usual  problem  of 
computing  the  convex  hull.  We  review  some  of  the  inclusion  and  exclusion  problems  that 
have  been  studied. 

(1)  We  mainly  focus  on  recent  results  within  the  milieu  of  computational  geometry. 
However,  it  should  be  pointed  out  that  there  is  a  related,  much  larger  and  older 
literature  arising  from  the  field  of  operations  research.  It  should  be  clear  that  our 
problems  arc  closely  related  to  the  "stock-cutting  problems"  which  arc  concerned 
with  cutting  a  sheet  of  material  into  smaller  subparts  under  various  constraints  (such 
as  all  subparts  are  congruent  to  a  given  shape)  and  subject  to  some  optimality 
criteria.  See  [DB83]  and  the  references  therein.  As  pointed  out  in  [DB83]  the 
enclosure  problems  is  a  key  subproblem  in  the  more  general  stock-cutting  problems. 
More  generally,  the  enclosure  and  inclusion  problems  can  be  viewed  as  polygon 
approximation  problems.  For  example,  for  the  purposes  of  detecting  collision  in 
robotics  we  typically  approximate  a  complicated  shape  by  a  simpler  enclosing  body. 


'  Alternatively,  these  might  be  called  the  inscription  and  the  circumscription  problems. 


(2)  The  potato-peeling  problem  is  the  case  Inc(Pau,  Peon'  area)  where  Pau  is  the  family 
of  all  simple  polygons  and  P^on  is  the  family  of  convex  polygons.  It  turns  out  that 
we  can  also  solve  the  potato-peeling  problem  in  the  case  where  perimeter  rather 
than  area  is  the  measure.  A  variation  of  the  potato-peeling  problem  docs  not  fall 
under  the  above  notion  of  inclusion  problems:  find  the  largest  convex  subregion  Q 
of  the  given  P  subject  to  the  constraint  that  Q  is  obtained  from  P  by  at  most  k  cuts. 

(3)  Dobkin  and  Snyder  [DS79]  considered  the  inclusion  problem  Inc{Pcon>  P^f  i^rea) 

I  where  Pj  denotes  the  class  of  all  triangles  and  in  general  /\  denotes  the  class  of  all 
convex  it-gons.  Their  algorithm  runs  in  linear  time.  This  result  was  extended  by 
Boyce,  et.  al.  [BDDG82]  to  the  problems  Inc(Pconf  Pkr  <^rea)  and 
lnc(Pcg„,  Pjc,  perimeter)  for  any  fixed  k.  The  running  time  of  these  algorithms  is 
0(k  n  log^n).  Note  that  unlike  the  potato-peeling  problem,  the  finiteness  of  these 
problems  is  easy  to  show:  it  follows  from  the  fact  that  the  vertices  of  any  maximal 
k-gon  must  be  a  subset  of  the  vertices  of  the  input  polygon.  The  techniques  of 
[DS79]  and  [BDDG82]  are  not  sufficient  for  the  more  general  problem  of 
Inc(Paii,  Pic,  area),  since  they  rely  on  the  convexity  of  the  input  polygon. 

(4)  Klee  and  Laskowski  [KL83]  considered  the  enclosure  problem  Enc(Pau,  Pj,  area) 
and  derived  an  Oin  log  /z)  solution.  O'Rourke,  et.  al.  [OAMB84]  improved  it  to 
linear  time.  De  Pano  [DP84]  described  how  the  method  in  [KL83]  extends  to  solve 
EnciPaii,  Pky  area)  for  all  ^  in  0(n  log  n)  time  (which  is  exponential  in  k). 
Chang  and  Yap  [CY84]  improved  De  Pano's  result  to  0(n^  log  k).  By  further 
refinement,  we  obtain  the  bound  of  0{n  log  n  log  k)  in  [ACY85].  We  remark 
that  [DB83]  claims  to  have  a  linear  time  solution  to  this  problem.  However,  their 
optimality  proof  is  faulty;  indeed  O'Rourke  [OR84b]  has  provided  some  counter 
examples. 

(5)  For  the  problem  of  finding  the  largest  rectangle  containing  a  given  polygon, 
Toussaint  [Tou79]  improved  a  previous  quadratic  time  solution  to  linear  time.    In 

general,  let  0  =  (9j,  .  .  .  ,  9j^  be  any  sequence  of  angles  with  each  9;  <  it  and 
k 

2  ^i—  ik—2)Ti.  Let  Pq  denote  the  family  of  convex  /t-gons  whose  interior 
1  =  1 

angles  are  given  by  the  sequence  0.  Thus  we  have  the  problem  of  finding  smallest 
polygon  from  Pq  enclosing  a  given  convex  polygon:  Enc(P^g„,  Pq,  p,).  O'Rourke 
and  his  co-workers  have  solved  several  of  these  problems.  For  the  3-dimensional 
versions  of  these  problems,  [OR84a]  describes  an  0(n  )  algorithm  for  the  smallest 
rectangular  box  enclosing  a  polyhedron. 

(6)  And  finally  problems  such  as  finding  the  smallest  square  containing  a  given 
polygon,  are  also  interesting.  In  general,  let  P  be  any  polygon  and  shape(P)  be 
the  family  of  polygons  obtained  by  the  transformations  of  scaling,  rotation, 
translation    and    reflection    of    P.     Then    we    have    the   fixed    shape    problems 


InciPaii,  shape(P),  p.)  and  EnciPau,  shape(P),  y.)  for  any  convex  P. 

2.   Preliminary 

For  the  rest  of  this  paper,  unless  otherwise  stated,  we  assume  a  fixed  but  arbitrary 
polygon  P  =  (vq,  Vj,  .  .  .  ,v„_i)  with  n  corners.  P  has  k  ^  0  reflex  corners: 
v,-^,  V;  ,  .  .  .  ,v,-^  (0  ^  /i  <  J2  <    ■  ■  ■    "^  ")•    Write  Uj  for  V/^. 


Vn-l 


Figure  1.  Polygon  /'  =  (vq,  .  .  .  ,v„_i) 


Given  P,  we  are  to  find  any  maximum  area  convex  subset  Q  contained  in  P;  it  is 
intuitively  clear,  but  rigorously  proved  in  [G008I]  ,  that  Q  is  a  convex  polygon.  For 
instance  if  P  is  convex  then  Q  is  unique  and  equal  to  P . 

First  we  introduce  some  notations.  A  chord  of  P  is  a  maximal  line  segment  fully 
contained  in  P  (note  that  there  could  be  line  segments  in  P  with  both  end  points  on  the 
boundary  of  P  which  are  not  chords).  A  chord  is  extremal  if  it  contains  two  or  more 
corners  of  P.  In  particular,  an  edge  of  P  is  always  contained  in  an  extremal  chord.  It  is 
clear  that  a  maximum  area  convex  polygon  must  be  the  intersection  of  P  and  k  half- 
planes  defined  by  k  chords  of  P:  Cj,  C2,  ■  ■  ■  .C^  which  pass  through  the  reflex  corners 
Ui,  U2,  .  .  .  .Ujt  respectively.  For  any  chord  C  of  P  passing  through  a  unique  reflex 
corner  u,  let  C  denote  the  closed  half-plane  determined  by  C  such  that  for  a  sufficiently 
small  disc  D  centered  at  u,  we  have  DHC^  C  P.    For  a  chord  C  that  passes  through 

more  than  one  reflex  corner,  the  context  will  make  it  clear  which  half-plane  is  intended. 

k 

Thus  the  convex  polygon  determined  by  the  chords  Cj,  C2,  ...  ,Ci^  is  PniClC^j. 

1  =  1 

Clearly  some  of  the  C,'s  may  be  redundant. 

We  first  answer  a  simple  case  of  the  problem  where  the  given  polygon  P  has  just 
one  reflex  corner  Wj.  From  the  above  observation,  the  problem  amounts  to  determining 
the  chord  C  through  u^  that  maximizes  the  area  of  PDC'*'.    Let  H  be  the  set  of  extremal 


chords  through  u^.    Clearly  |H|    <  n.    Consider  the  butterfly  region  B  determined  by  a 
pair  of  adjacent  chords  au^a'  and  buib'  taken  from  H: 


Figure  2.  A -butterfly  with  center  Uj 

We  call  Uj  the  center  of  the  butterfly  and  the  line  segments  [a,  b]  and  [a' ,  b'\  the 
tips  of  the  butterfly.  A  chord  C  embedded  in  B  is  determined  by  any  point  c  in  the  tip 
[a,  b]  together  with  c'  in  [a' ,  i']  such  that  c,  Uj  and  c'  are  colinear;  C  =  cu^c' .  The 
problem  reduces  to  choosing  for  each  butterfly  a  chord  C  embedded  in  it  such  that  the 
area  a  =  C  PlS  is  maximized,  a  is  either  the  union  of  triangles  Awj  be  and  Auj  a'c' 
or  the  imion  of  the  triangles  Auj  ac  and  Au^  b'c' ,  depending  on  the  orientation  of  C. 
The  following  lemma  shows  that  either  C  is  extremal,  or  Uj  is  the  midpoint  of  the  chord: 

l"ic|  =  |"ic'|. 
In  the  later  case,  we  have  a  balanced  chord. 


\  / 


Figure  3.  A-  and  V-Butterflies 


The  pair  (Lg;  Lj)  of  the  lines  through  the  two  tips  of  B  are  the  supporting  lines  of 
B.  If  Lq  and  Lj  intersect  at  a  point  o  and  o  ^  C*"  (pick  any  C  embedded  in  B,  say 
[a,  a']),  then  we  call  the  butterfly  an  A -butterfly.  Otherwise  it  is  a  V -butterfly.  Sec 
Figure  3.    (See  also  the  next  section  for  the  general  context  for  the  'A'  and  'V  notations.) 


Lemma  1.  (Butterfly  Lemma)  Given  the  butterfly  B  determined  by  an  adjacent  pair 
of  extremal  chords  au-^  a'  and  bu\  b' ,  let  C  =  c«i  c'  be  a  chord  embedded  in  B 
maximizing  a.  If  fl  is  an  A-butterfly  then  C  is  either  balanced  or  extremal.  Otherwise 
S  is  a  V-butterfly  and  C  is  extremal. 

Proof.  Consider  the  case  where  B  is  an  A-butterfly.  First  without  loss  of 
generality,  assume  \aui\  ^  \ci'ui\.  If  \bui\  2:  l^'u^l,  then  by  a  simple  continuity 
argument  there  is  a  unique  balanced  chord  C*  embedded  in  5.  It  is  easy  to  see  that  any 
other  chord  C  determines  a  smaller  area  than  C*  (see  also  [CL71]).  Otherwise 
\bui\  <  \b'ui\  and  it  is  not  hard  to  see  that  the  extremal  chord  [b,  b']  maximizes  a.  In 
the  case  of  a  V-butterfly,  if  C  is  not  extremal  we  can  perturb  C  so  as  to  enlarge  a. 
Hence  we  conclude  that  C  must  be  equal  to  one  of  the  two  extremal  chords,   n 

This  lemma  clearly  leads  to  a  linear  time  algorithm  for  the  potato-peeling  problem  if 
P  has  one  reflex  corner.  In  the  next  section,  we  look  at  the  general  case  and  at  A-  and 
V-butterflies  in  a  more  general  setting. 

3.    Series  of  Butterflies  and  Chains  of  Chords 

In  this  section  we  will  give  a  finiteness  criterion  for  the  potato-peeling  problem. 
First  let  us  introduce  the  terminology  related  to  "sequences  of  butterflies"  and  "sequences 
of  chords."  It  turns  out  that  sequences  of  butterflies  can  be  classified  into  two  types  with 
quite  different  algorithmic  properties.  For  simplicity,  we  assume  that  no  three  corners 
are  colinear.  As  in  the  previous  section,  for  each  reflex  corner  u  of  the  polygon  P,  we 
can  form  a  circular  list  of  all  extremal  chords  through  u.  A  butterfly  of  P  is  the  region 
determined  by  a  pair  of  adjacent  chords  through  u  having  the  obvious  shape.  Our 
definitions  for  butterflies  and  chords  are  relative  to  some  fixed  polygon  P. 

Definition.  A  series  of  butterflies  S  =  (Sj,  .  .  .  ,B„),  m  ^  1,  is  any  sequence 
of  butterflies  satisfying 

(1)  If   c,-   is    the    center    of   S,-   then    (cj,  C2,  ■  ■  ■  ,c„)    forms    a   convex    polygon    Q 
contained  in  P.    Q  is  degenerate  if  m  =  1. 

(2)  The   two   wings   of  each   butterfly    are   ordered   so   that  the   'forward'   wing  of  5j 
intersects  the  'backward'  wing  of  5,  +  ^  for  /  =  1,  .  .  .  ,/n  —  1. 


Let  Cq  (resp.  Cj)  be  the  tip  of  the  backward  (rcsp.  forward)  wing  of  5^  (resp. 
Bn).  Then  (Cq,  C {)  is  called  the  (pair  of)  supporting  tips  of  the  series.  If  L,-  is  the 
line  through  C (  then  (Lq,  L{)  is  the  pair  of  supporting  lines.  Note  that  it  is  possible  for 
Lq  and  Li  to  be  parallel  or  even  equal.  Let  m  >  1.  If  the  pair  of  supporting  lines  are 
coincident  or  parallel,  or  if  they  intersect  at  a  point  on  the  side  of  the  line  c^c^  opposite 
to  Q,  we  will  say  (Lq,  L{)  and  (Cq,  C {)  are  V-shaped.  Otherwise  we  say  they  are  A- 
shaped.  If  a  series  is  supported  by  a  pair  of  V-shaped  lines,  then  it  is  a  V-series. 
Otherwise  it  is  an  A-series .  These  definitions  are  seen  to  be  extensions  of  the  previous 
definition  of  an  A-  or  V-butterfly.  (Remark:  'V  and  'A'  are  chosen  for  the  shapes  of 
these  letters,  being  mnemonic  for  the  orientations  of  the  supporting  lines  -  if  we  imagine 
the  line  CiC„  as  horizontal  and  Q  as  sitting  above  c^c^.) 


Figure  4.  A  V-chain  and  an  A -chain 


Notation:  Let/?  \q  denotes  the  reflection  oi  p  about  q,  i.e.  p  |^  is  a  point  on  the  line 
pq  such  that  q  is  the  midpoint  between  p  and  p  \q. 

Now  we  present  some  definitions  for  chords.  Let  (Cg,  C^+i)  be  a  pair  of  chords, 
not  necessarily  distinct.  A  (Cq,  C„  +  i)-chain  of  chords  is  a  sequence  of  chords 
(Cj,  C2.  ■  ■  ■  .C„)  such  that  C,  intersects  C,+  i  at  x,-  for  i  =  Q,  .  .  .  ,m  and 
Q  =  (xq,  Xi,  .  .  .  ,x„)  is  a  convex  polygon.  If  (Cq,  C„  +  {)  is  understood,  we  just  say 
'chain'.    Q  is  called  the  core  of  the  chain.    Call  [x,_i,  xj  the  truncated  version  of  the 


chord  Ci  (f  =  1,  .  .  .  ,m).  The  x/s  are  the  nodes  of  the  chain.  The  chain  is  said  to  be 
balanced  if  for  each  /  =  \,  .  .  .  ,m  the  truncated  version  of  C^  passes  through  reflex 
corners  p,-,  qi  (possibly/)/  =  Qi)  and  has  the  following  properties: 

(1)  pi  and  qi  both  lie  in  the  truncated  chord  [j:,_i,  xj.  We  may  assume  that  /?,•  lies 
between  x,_i  and  ^,-. 

(2)  The  midpoint  of  [j:,_i,  xj  lies  in  [p,,  <jj.    So  if  p,  =  ^,-  then  p,-  is  the  midpoint. 

(3)  (^Bracketing  property)  x,!^,-  lies  between  jc,_i  and  J:,-ilp/.  This  property  is  so- 
called  because  we  imagine  x^lqi  to  be  a  left  bracket  and  x,_i|p,-  to  be  a  right 
bracket.    Note  that  if  p,-  =  <y,-  then  J:,_i|pi  =  J:,-  and  x,|^,-  =  Xi-i- 


Pi 


•^.ki     ^,-llPi 


Figure  5. 


If  pi  ¥=  qi  we  call  (p,-,  ^j)  a  double-pivot,  otherwise  it  is  a  single-pivot.  If  the 
balanced  chain  has  only  single-pivots  then  it  is  a  simply  balanced  chain.  If 
S  =  (5i,  B2,  ■  ■  ■  >B„)  is  a  series  of  butterflies,  then  a  chain  embedded  in  B  is  simply 
a  (Co,  C„+i)-chain  C  =  (Ci,  C2,  .  .  .  XJ  where  (Cq,C„  +  i)  is  the  pair  of 
supporting  tips  of  B  and  each  C,-  is  embedded  in  fi,-.  Note:  If  C;  has  a  single-pivot  p,- 
then  p,  is  the  center  of  5,;  otherwise  C,  is  necessarily  one  of  the  two  extremal  chords  of 
5,-.    The  concept  of  double-pivots  is  not  relevant  until  the  next  section. 

Let  B  =  (Bi,  B2,  ■  ■  ■  ,B^  be  a  series  and  (Cj,  C2,  .  .  -  .C „)  be  a  sequence  of 
chords    such    that    each    C,    is    embedded    in    S,.      Here    we    do    not    assume    that 

(Cj,  C2 C „)   is  necessarily  a  chain  (i.e.  some  C,-,  C,  +  i  may  not  intersect).    We 

say  the  sequence  is  optimal  (for  the  series  it  is  embedded  in)  if  the  area  of 

m 

pn(n  ct) 
1  =  1 

is  a  local  maximal,  i.e.,  any  sufficiently  small  perturbation  of  the  chords  produces  a 
sequence  with  smaller  area.  A  sequence  of  chords  is  optimum  if  the  area  is  maximum 
over  all  chains  embedded  in  the  series  of  butterflies.  If  the  optimal  or  optimum  sequence 
of  chords  turns  out  to  be  a  chain,  we  call  it  an  optimal  or  optimum  chain. 


10 


Before  proving  a  basic  lemma  next,  we  introduce  a  useful  concept.  Let  Lq  be  a 
fixed  line  and  let  V  be  the  corners  of  polygon  P  and  {  pi,  .  .  .  ,p„]  he,  some  subset  of 
the  reflex  corners  of  P.   An  (Lq,  pj,  .  .  .  ,p^-path  is  a  polygonal  path  of  the  form 

n  =  (xq,  Xi xj 

where  Xq  €  Lq  and  for  /  =  1,  .  .  .  ,m  and  p,-  is  the  midpoint  of  [j:,_i,  a:,].    For  any 
point  X  and  line  L,  let  L(x)  denote  the  line  through  x  parallel  to  L. 


LQiX2) 


X2 


Lo(Xi) 


Figure  6. 


Note  that  the  line  Lq(x^  is  a  function  of  Lq  and  P\,P2,  ■  ■  .  ,Pi  only  (i.e.  Lq{x^  is 
independent  of  x,  or  of  the  particular  choice  of  IT).  Furthermore  for  any  z  =  0,  .  .  .  ,m, 
any  point  x  on  Lq{xI)  determines  a  unique  {Lq,Pi,  .  .  .  ,/?^-path  and  vice-versa.  In 
particular,  the  choice  of  any  x  i  Lq  determines  a  path.  We  can  think  of  a  path  as  a 
configuration  in  a  system  of  "interconnecting  levers;"  each  [x,_j,  xj  is  an  stretchable 
lever  on  the  fixed  axis  pi  and  nodes  x,_i  and  x,-  are  constrained  to  glide  along  parallel 
slots  Lo(x,_i)  and  LqCx;)  respectively.  We  call  Xq  ^  Lq  critical  if  the 
(Lq,  p  I,  .  .  .  ,p^)-path  at  Xq  has  the  property  that  for  some  i, 

(a)  either  [x,_i,  xj  passes  through  some  corner  in  V—{pi,  .  .  .  ,p„  } 

(b)  or  x,_i,  X,-,  x,>i  are  colinear. 

A  minimal  interval  I  Q  Lq  bounded  by  two  critical  points  is  called  a  critical  interval. 
An  {Lq,  Pi,  .  .  .  ,p^-path  is  said  to  belong  to  /  if  its  first  node  Xq  is  in  /.  If  the  nodes 
(-'^O'  -"^l-  •  •  •  •^m)  forms  the  corners  of  a  simple  polygon  Q  in  the  indicated  order  then 
we  define  the  area  of  the  path  to  be  the  area  of  Q.  Then  the  area  of  Q  depends  only  on 
(^0.  P  l>  ■  •  •  •Pm)  ^^^  °o'  o°  n.    This  area-invariance  property  is  due  to  the  fact  that 
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each  Pi  is  the  midpoint  of  the  segment  [jc<_i,  xj,  /  =  0.  .  ,  .  ,m  and  hence  as  the 
"levers"  go  up  and  down,  Q  loses  exactly  as  much  area  as  it  gains.  We  exploit  this 
property  to  show 

Lemma  2.  (V-Lemma)  Let  C  =  {C i,  C2,  ■  .  .  ,Cm)  ^^  *  chain  that  is  optimal  for  a 
V-scries  of  butterflies  B  =  (Sj,  Bj,  ■  ■  ■  .B^.  If  the  supporting  lines  of  B  are  not 
parallel  or  coincident  then  at  least  one  of  the  C,'s  is  extremal.  If  the  supporting  lines  are 
parallel  or  coincident  then  we  can  modify  the  chords  without  decreasing  the  area  defined 
by  the  chain  so  that  at  least  one  of  the  chords  becomes  extremal. 

Proof.  The  case  m  =  1  is  the  consequence  of  the  Butterfly  Lemma.  So  let  m  >  1. 
For  each  1  S  i  ^  m,  consider  the  butterfly  B  centered  at/j^  and  supported  by  C,_i  and 
C,-.  According  to  Butterfly  Lemma,  Cj  is  eithei  extremal  or  balanced  when  truncated.  If 
Ci  is  extremal  then  we  are  done.  Thus  assume  without  loss  of  generality  that  the 
truncated  version  of  each  chord  C/  (/  =  1,  .  .  .  ,m)  is  balanced.  However  if  Lq  and  L^ 
are  not  parallel,  we  can  perturb  the  chord  C „  out  of  balance  while  keeping  all  other 
chords  balanced,  resulting  in  another  cLaiu  C'  ~  {Ci,  C2,  •  ■  •  ,C„)  with  nodes 
(xq,  xi,  .  .  .  .x„).  (See  Figure  7.) 


^0(-^2) 


Figure  7.  The  dotted  chain  represents  the  perturbation 


Hence  (xq,  Xj J^m-l)    is  an  (Lq.Pi.Pi-  •  •  •  .Pm-l)-path  and  C„  and  C„  are 

embedded  in  the  V-shaped  butterfly  supported  by  Lq(x„)  and  Li(x„  +  i)  =  Lj.  li  m  = 
even,  and  the  perturbation  Xq  is  in  a  direction  along  Lg  away  from  the  intersection  of  Lg 
and  Li  then  the  area  of  C  is  increased.  If  m  =  odd,  we  perturb  in  the  opposite 
direction.  More  precisely,  in  addition  to  the  area-invariancc  of  the 
(Lq,  p  I,  p2.  ■  ■  ■  ,P;„_i)-path,  C„  contributes  a  gain  in  area  that  equals  the  area  of 
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A  X2X2y  in  Figure  7.  Thus  C  results  in  a  larger  area  than  C,  contradicting  the  local 
optimality  of  C-  Finally  if  the  lines  are  parallel  or  coincident,  we  can  perturb  x„  along 
Lq(x„)  =  Li  without  changing  the  area  of  the  series  since  the  perturbed  series  is  an 
(Lq,Pi,  .  .  .  ,p„)-path.  We  can  perturb  x„  until  the  corresponding  Xq  is  critical, 
implying  that  one  of  the  chords  is  extremal.   D 

We  next  investigate  the  considerably  more  subtle  A -series  of  butterflies. 

Lemma  3.    (A-Lemma)  Let  S  be  an  A-series  of  butterflies.  Then 

(a)  B  has  at  most  one  simply  balanced  chain  C*  embedded  in  it. 

(b)  An  optimal  chain  of  B  is  the  unique  simply  balanced  chain  C*  when  C*  is  embedded 
in  B,  otherwise  the  chain  has  an  extremal  chord. 

Proof,  (a)  is  proved  in  the  next  section  in  a  slightly  more  general  setting.  To  see 
(b),  suppose  C,-  is  an  unbalanced  chord  in  a  chain  with  no  extremal  chords.  The  same 
argument  in  V-Lcmma  applies  here:  perturbing  Cj  toward  the  balance  position  (this  is 
possible  since  C,-  is  not  extremal),  while  keeping  all  other  chords  unchanged,  increases 
the  area  of  the  series  and  therefore  violates  the  local  optimality.   a 


Lemma  4.  Let  Q  be  a  maximal  convex  polygon  in  P  bounded  by  the  sequence 
(Cq,  C I,  .  .  .  .C„)  of  chords.  Then  Q  contains  at  least  two  extremal  chords  of  P. 
Furthermore,  there  exists  a  maximal  convex  polygon  with  three  or  more  extremal  chords. 

Proof.  To  show  the  first  part,  assume  to  the  contrary  that  Q  has  0  or  1  extremal 
chord.  Note  that  (C2,  .  .  .  ,C„)  forms  a  (Cg,  Ci)-chain.  Without  loss  of  generality, 
assume  that  this  chain  contains  no  extremal  chords.  Clearly  (C2,  .  .  .  ,C„)  is  a  V-chain, 
and  by  the  V-Lemma,  it  is  not  optimal.  For  the  second  part,  note  that  if  there  are  only' 
two  extremal  chords,  say  Cq  and  Cj^,  then  they  must  be  parallel  (otherwise  the  V-Lemma 
shows  that  either  the  (Cq,  Cjt)-chain  or  the  (C^,  Co)-chain  is  non-optimal).  But  in  that 
case,  the  proof  of  the  V-Lemma  shows  how  to  perturb  the  (Cq,  Cj) -chain  to  get 
additional  extremal  edges  while  preserving  maximality.    D 


The  above  lemmas  provide  us  with  a  finiteness  criterion  in  the  sense  that  we  can 
guess  that  Q  is  determined  by  a  sequence  Cj,  C2.  ■  ■  ■  ,Cj^  of  extremal  chords  together 
with  series  of  butterflies  supported  by  (C,-,  C,+  i)  for  /  =  1,  .  .  .  ,yt.  This  gives  an 
exponential  time  algorithm  provided  that  we  can  find  simply  balanced  chains  for  any 
given  series  of  butterflies  in  exponential  time.  We  next  show  that  such  chains  can  in  fact 
be  found  in  polynomial  time. 
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4.    A  Geometric  Problem 

The  problem  of  finding  the  balanced  chain  in  an  A -series  of  butterflies  can  be 
reduced  to  an  abstract  geometric  problem.  First  we  transplant  some  notations  from  the 
previous  section  to  a  different  geometric  setting.  We  assume  a  fixed  set  R  of  points  in 
the  plane.  Let  L"*"  be  the  half-plane  to  the  right  of  a  directed  line  L.  Given  a  pair  of 
directed  lines  (Lq,  Li),  let  a  chord  denote  a  line  segment  contained  in  Lq  flLj  passing 
through  at  least  one  point  of  R  and  with  endpoints  in  Lq  and  Lj.  A  chord  is  extremal  if 
it  passes  through  two  or  more  points  of  ^.  As  before  the  points  of  /?  in  a  chord  arc 
called  the  pivots  of  the  chord.  The  definitions  of  chains  and  nodes  are  the  same  as  in 
Section  3.  Let  C  =  (Ci,  Ci,  .  .  .  ,CJ  be  an  (Lq,  Li)-chain  with  {xq,  x^,  .  .  .  .x^  as 
nodes. 


Figure  8.  (Cj,  Ci,  Cy)  is  an  empty  (Lq,  Lj) -chain  (the  points  of  R  are  indicated  by  'x  ') 

The  convex  polygon  Q  =  {xq,  x^,  .  .  .  ,X;„)  is  the  core  of  C-  The  polygonal  path 
{xq,  Xj,  .  .  .  ,x^)  partitions  Lq  r\Li  into  a  finite  and  an  infinite  regions.  Let  p(C) 
denote  the  open  infinite  region  so  defined.  We  say  C  is  empty  with  respect  to  R  if  p(0 
does  not  intersect  R.  The  definition  of  balanced  chains  is  the  same  as  in  the  previous 
section,  with  the  set  R  playing  the  role  of  the  corners  of  P.  (As  usual;  assume  for 
simplicity  that  no  three  points  of  R  are  colinear.)  Now  let  us  consider  the  following 
problem: 

Let  R  be  the  given  set  of  points  on  the  plane.  For  each  pair  of  directed  lines 
{Lq,  Li)  such  that  Lq  and  Lj  each  passes  through  an  ordered  pair  of  points  from  R, 
find  the  balanced  (Lq,  Lj) -chain  that  is  empty. 


Theorem  5.    Let  (Lq,  Li)  be  a  pair  of  directed  lines  and  R  he  a  set  of  points  on  the 
plane.    Then  there  is  a  unique  balanced  (Lq,  Z,i)-chain  that  is  empty  with  respect  to  R. 

We  show   a  useful  lemma   along  the  way  to  proving  this  theorem.     The   following 
notations  are  needed  for  this  lemma: 
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Notations:  Let  C  be  a  balanced  chain  and  let  (p  i,  q^  be  the  pivot  of  the  chord  C i  in 
C  (possibly  Pi  =  qd  and  let  (xq,  Xi,  .  .  .  ,x^  denote  the  nodes  of  C.  li  p i  ¥=  q^  then 
assume  ^,-  lies  between  p,-  and  Xi-  Recall  that  for  any  point  p,  L\{p)  denotes  the  line 
through  p  and  parallel  to  Lj.  Define  the  region  Si  to  be  the  strip  between  the  parallel 
lines  Li(j:,_i)  and  L\(^Xi).  Define  the  region  A,  to  be  the  strip  between  the  parallel  lines 
Li(j:,_i)  and  L^ipi):  it  is  important  to  note  that  we  exclude  the  line  Li(j:,_i)  from  A,- 
but  include  Li(pi)  in  A;.  The  region  5,-  is  the  strip  between  L^ipi)  and  Li(x,_i|p ;): 
again  L^ip^  is  excluded  but  Li(jc,_i|pj)  is  included.  Similarly,  D,-  (resp.  E,)  is  the  strip 
and  between  the  lines  Li(x,|^j)  and  Liiqi)  (resp.  Li(q^  and  Liix^).  Recall  that  the 
chain  C  divides  the  quadrant  Lq  HLi  into  a  finite  region  and  an  infinite  p(C).  It  is 
convenient  to  regard  this  finite  region  to  be  above  the  chain  and  p(C)  to  be  below.  Thus 
each  of  the  strips  5,-,  A,-,  5,-,  etc.,  will  be  broken  up  into  two  half-strips,  above  and 
below  the  chain,  respectively.    See  Figures  9  and  10  for  these  regions. 


Li(x,_i)         LiipD  Li(Xj) 


Xi-l 


Pi^li 


Ai=Di     I     Bi=Ei 


Xi-l 


Di  ,    Ei 


U- 


I      I 
I      I 


A,     .     Bi 


'  ^i 


Figure  9.  The  regions  A,,  D,,  Bi  and  E,-  (two  cases) 

Let  a,-  =  A,UD,-  and  P,-  =  BiUEi-    Note  that  the  a,-  and  P,'  are  non-disjoint  because  of 
the  bracketing  property  and  thatp,-  =  q^  iff  A,-  =  D,  iff  5,-  =  E,-. 


Lemma  6.  Let  Lq,  Lq,  L^  be  distinct  lines.  Let  C  =  (Cj,  C2.  ■  ■  ■  ,C^  (resp. 
C  =  (Cj,  C2,  ■  ■  ■  ,C„'))  be  a  balanced  {Lq,  Li)-chain  (resp.  (Lq,  Li)-chain).  Let 
Xj'%,  pj's  be  the  nodes  and  pivots  of  C  .  If  both  C  and  C  are  empty,  then  Xj  cannot  lie 
above  C  in  a,-  for  any  i  and  j. 
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Figure  10.  The  shaded  areas  are  forbidden  for  nodes  of  C:  the  two  cases 


Proof.  When  j  =  m' ,  by  definition  x„'  lies  in  E,  and  hence  is  not  in  any  a,. 
Choose  j  to  be  the  largest  index  such  that  Xj-i  (j  ^  m')  is  above  the  chain  in  a,-  for 
some  I.  Hence  Xj  is  lying  below  C  or  not  in  o,-.  Wi  will  prove  that  Xj  or  Xj  +  i  lie  above 
the  chain  in  some  a^  (k  >  i),  thus  contradicting  our  choice  of  y.  Surely  Xj  is  not  above 
the  chain  in  3i  for  otherwise  p,-  would  be  in  p(C')  contradicting  the  claim  that  C  is 
empty.    So  Xj  is  either  below  the  chain  in  5,-  or  in  5^  for  some  k  >  i. 

First  consider  the  case  where  X:  is  in  5^.  If  Xj  lies  below  the  chain,  then  pj  would 
be  below  C-    This  contradicts  the  assumption  that  C  is  empty. 
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Li(x^ 


Figure  11. 


Suppose  Xj  lies  above  the  chciin  in  E^— a^  (the  proof  for  5^—0^  is  similar).  Cj  must 
intersect  L\{qj^  below  the  chain  (otherwise  q^^  would  be  in  p(C').  contradicting  the 
emptyness  of  C).  Let  z  denote  this  intersection  point.  Similarly  C,  intersects  L\{pi) 
below  C I  at  some  point  y  (see  Figure  11).  Then  the  midpoint  of  the  segment  [x,_i,  Xj] 
lies  in  \y ,  z].  Note  that  p.-  must  lie  to  the  right  of  y  and  qj  to  the  left  of  z  to  satisfy  the 
bracketing  property.  Since  p,  must  lie  to  the  left  of  q;,  this  implies  both  p,  and  q:  lie 
below  the  chain  C,  contradiction.  Thus  x,  must  lie  in  some  a^.  This  contradicts  our 
choice  of  the  index  j.  Finally,  consider  the  case  when  X:  stays  in  5,-  below  the  chain  C. 
Similar  arguments  apply  to  show  that  -Cj  +  j  is  above  the  chain  in  aj^.  This  contradiction 
completes  the  proof,   n 
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Lq        Li{Xq) 


^l(Pl)  ^lUolPi) 


Figure  12. 


Proof  of  the  theorem:  There  are  two  parts  to  this  theorem:  that  there  cannot  be  more  than 
one  balanced  chain  and  that  there  exists  at  least  one  balanced  chain. 

(Uniqueness)  Let  C  =  (Cj,  C2,  ■  ■  ■  .C,,.)  be  a  balanced  chain.  We  derive  a 
contradiction  by  assuming  the  existence  of  another  balanced  chain 
C  =  (Ci,  C2.  ■  ■  •  .C„').  Let  pi'&,  X:'i  and  /?,'s,  Xj's  be  the  pivots  and  nodes  of  C 
and  C  respectively.  To  apply  the  previous  lemma  we  only  have  to  show  that  x^  lies  in 
aj  for  some  i,  j.  Initially  assume  Xq  #  Xq.  Suppose  Xq  is  below  xq  on  Lq.  (See  Figure 
12.)  Observe  that  JCj  must  lie  above  C  otherwise  Pi  is  below  C.  li  Xj  lies  in  any  a,-  we 
are  done.  So  let  x\  lie  in  (3,- a,-.  Note  that  x[  is  to  the  right  of  LiixQ\p{).  li  Xi  is  in 
5^  then  p[  is  below  C,  contradiction.  If  x'l  is  in  Ej  — Dj,  let  z  be  the  intersection  of  Cj 
with  LqC^i).  Then  z  must  be  below  C.  Note  that  the  midpoint  of  Cj  lies  to  left  of  z  and 
hence  pi  lies  to  the  left  of  z.  This  implies  pj  is  below  C,  contradiction.  Therefore  Xq 
cannot  lie  below  Xq.  If  xq  is  above  Xq  on  Lq  then  it  is  easy  to  show  that  Xq  is  in  Di 
(otherwise  cither  p  j  is  in  p(C')  or  C  never  reaches  Lj). 

It  remains  to  consider  the  possibility  Xq  =  Xq.  Since  C  ^  C  ,  let  Xj  be  the  first 
node  such  that  jc,  ^  JC..  Suppose  Xj  does  not  lie  in  the  line  through  [x^_i,  Xj].  An 
analysis  similar  to  the  above  shows  that  Xj  lies  in  ajt  for  some  k.  Finally  suppose  Xj  lies 
in  the  line  through  [xj^i,  Xj]:  if  it  lies  in  the  segment  [.ry_i,  Xj]  then  Xj+i  must  lie  in 
some  Ojt  by  the  same  argument  as  the  previous  case.  Otherwise,  Xj  clearly  is  in  Aj+i  and 
again  we  have  a  contradiction. 

(Existence)  The  existence  proof  can  be  regarded  as  an  algorithm  although  we  do  not 
know  of  a  polynomial  time  bound  on  its  complexity. 
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Wc  now  give  a  'scan  line'  algorithm  for  computing  the  balanced  (Lq,  Lj) -chain.  A 
sequence  of  (Lo(^),  Li)-chain8  are  computed  where  LqCO  is  a  line  parallel  to  Lq  at  a 
distance  t  to  the  left  of  Lq.  Imagine  the  line  LqCO  moving  from  the  infinite  left  toward 
Lq  as  t  approach  0  from  oo.  As  r  -»  0,  LqCi)  -  Lq  and  the  (LqCOi  Li)-chain  becomes  the 
(Lq,  Li)-chain.  During  the  process,  there  are  events,  Iq,  ti,  .  .  .  .t„  that  divide  the 
scanning  process  into  intervals  where  changes  in  the  (Lq^I) ,  Li)-chain  are  smooth  within 
each  interval.  More  specifically,  between  events  the  pivots  remain  the  same  and  the 
slope  of  each  chord  changes  at  a  smooth  rate  (with  respect  to  t).  Initially  the 
(Lo(°°),  Li)-chain  consists  of  an  infinite  chord  parallel  to  L^  through  the  point  pj  in 
RHLi  that  is  farthest  from  Lj.  Then  this  chord  turns  continuously  counterclockwise 
about  pi  as  LqC^)  gradually  moves  toward  Lq  until  the  chord  hits  a  new  point  qiin  R  and 
turns  into  a  double  pivot.    This  is  the  first  critical  moment  Iq. 

As  Lo(0  continues  to  move  closer  to  Lq  and  xq  moves  toward  Xj  (as  usual,  jc,- 
denotes  the  i'  node  of  the  balanced  chain,  C,  is  the  i'  chord,  etc.),  Xq\pi  moves  toward 
both  Xi\qi  and  pi  and  the  midpoint  rj  of  Cj  moves  toward  ^j.  Two  things  can  happen 
at  the  next  critical  moment  t  =  t-^:  either  (a)  Xq\pi  meets  x-y\qi  and  the  double-pivot 
chord  'splits'  into  two  single-pivot  chords  balanced  at  p^  and  q^,  or  else  (b)  rj  meets  q^ 
and  the  chord  starts  turning  counterclockwise  around  q^  'leaving'  p\  behind,  (a)  If  the 
splitting  event  occurs  first,  wc  will  have  a  chain  of  size  2,  and  as  LQ(t)  moves  further  in, 
in  order  to  balance  C^  and  Ci,  Xq  and  Xj  move  inwards  and  jTj  moves  outwards  along  the 
direction  of  Lj.  (Note:  'inwards'  and  'outwards'  are  with  respect  to  the  center  of  Q  of 
the  chain.)  This  continues  until  t  =  t2  when  a  point  in  R  is  hit  by  either  Cj  or  C2  and 
turns  a  single-pivot  chord  into  a  double-pivot  chord,  (b)  If  the  leaving  event  occurs  first 
the  double-pivot  chord  turns  into  a  single-pivot  chord  which  keeps  turning 
counterclockwise  until  f  =  ^2  it  hits  another  point  and  turns  itself  into  a  double-pivot 
chord.  So  this  process  of  forming  double  pivots,  splitting  and  leaving  continues  until 
t  =  t„  when  Lo(0  reaches  Lq. 

In  general,  the  (Lo(^),  Z-i)-chain,  (Cj,  C2,  ■  .  ■  ,C^  has  an  initial  sequence 
(possibly  empty)  of  single-pivot  chords  Cj,  .  .  .  ,Cj-i  and  the  leftmost  double-pivot 
chord  Cj  (if  there  are  no  double-pivots,  take  _/'  =  m  +  1  and  C,-  =  Lj)  followed  by  some 
other  chords.  At  r  =  t/^,  Lq^I)  moves  inwards  the  even-numbered  nodes  {X2i,  2/  <  j) 
move  inwards  while  odd-numbered  nodes  {x2i+\,  2z  +  1  <  j)  move  outwards  along  the 
direction  of  Cj.  Observe  that  C,  turns  clockwise  if  /  is  even  and  counterclockwise 
otherwise,  for  z  =  1,  .  .  .  ,j—\.  The  length  of  the  truncated  version  of  C,  increases  or 
decreases  depending  on  the  parity  of  y.  The  rest  of  the  (Lo(/') .  ^i)-chain  remains 
unchanged.  We  can  classify  the  events  that  could  happen  next  at  f  =  fj^  +  j  into  five 
categories: 
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(1)  (flattening)  Two  consecutive  chords  Cii  and  C2/+i  could  flatten  out  and  become  a 
double-pivot  chord.    It  should  be  noted  that  C2,_i  and  Cit  cannot  flatten  out. 

(2)  (hitting)  One  of  the  single-pivot  chord  could  hit  a  point  in  R  thus  becomes  a 
double-pivot  chord. 

(3)  (splitting)  If  j  is  odd,  the  truncated  version  of  Cj  could  shorten  to  such  an  extend 
that  ;c,_i|py  meets  Xj\qj.  Then  Cj  loses  the  bracketing  property  and  splits  into  two 
single-pivot  chords. 

(4)  (leaving  p.)  When  C,  is  shortening,  another  situation  could  also  arise.  The 
midpoint  r,-  could  meet  qj  before  splitting  occurs.  Subsequently,  Cj  will  turn 
counterclockwise  around  qj  leaving  Pj  behind.  In  effect,  Cj  turns  into  a  single- 
pivot  chord  balanced  at  qj. 

(5)  (leaving  qj)  Finally,  if;  is  even  the  chord  Cj  is  lengthening  and  the  only  event  that 
can  happen  is  the  midpoint  r,  moving  left  and  meeting  pj.  Subsequently,  Cj  will 
pivot  clockwise  around  pj  leaving  qj  behind 

The  algorithm  first  establishes  a  balanced  chain  for  (Lito),  L^)  and  then  repeats  the 
following  step  until  Lo(0  reaches  Lq: 

Find  the  leftmost  double-pivot  chord,  Cj  on  the  current  (Lo(^),  Li)-chain.  In  the  case  of 
a  chain  consisting  of  just  single-pivot  chords,  take  Lj  for  Cj.  It  is  not  too  hard  to 
compute  the  values  of  t  when  flattening,  leaving  or  splitting  occur  at  each  C,-,  since  those 
are  determined  by  the  chain  alone.  The  hitting  event  at  each  Cj  is  more  difficult  to 
compute  since  it  involves  points  not  on  the  chain.  But  it  is  clear  the  point  hit  by  C,-  has  to 
be  on  the  convex  hull  of  some  subset  of  R  containing  p,.  So  we  partition  the  plane  into; 
strips  and  two  half-planes  by  the  lines  Lj{t^,  i  =  0,  .  .  .  J  where  Ly(fj)  is  the  line 
through  Xf  parallel  to  Cj.  Let  //,•  be  the  convex  hull  of  those  points  of  R  in  the  strip 
between  Lj(ti-i)  and  Lj(t^.  The  point  hit  by  an  even-numbered  chord  C2,-  is  the  point 
on  H2i  clockwise  from  p2i;  the  point  hit  by  C2,  +  i  is  the  next  point  on  //2i  +  l 
counterclockwise  from  P2i+1-  ^°  ^°^  each  i  we  can  compute  the  value  of  t  when  an  event 
involving  C,-  occurs.  The  next  event  at  f  =  Z^  +  j  is  determined  by  one  with  the  smallest 
of  such  t  values.  To  complete  the  present  step,  we  just  update  (Lo(0»  ^i)-chain 
accordingly. 

This  concludes  our  proof  of  the  existence  of  a  unique  balanced  chain  for  any 
(Lq,  Li). 

S.    Decomposition  of  Balanced  Chains 

Let  /?  be  a  fixed  set  of  n  points,  and  {pj,  .  .  .  ,p„  }  Q  R-  Recall  the  definition  of 
(L.pi,  .  .  .  ,p „)-paths,  critical  points  and  critical  intervals  in  Section  3.  In  this 
section  we  use  these  concepts  to  describe  the  composition  and  decomposition  processes  of 
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balanced  chains. 

Lemma  7.  Given  (Lq,  p\,  .  .  .  ,p„^,  we  can  determine  all  the  critical  intervals  of 
(Lq,  Pi,  .  .  .  ,/7  J-paths  in  0(n  log  n)  time. 

Proof.  Pick  an  arbitrary  point  Xq  on  Lq.  Form  the  unique  path 
n  =  {xq,  x^,  .  .  .  ,x^  and  divide  the  plane  into  strips  determined  by  the  parallel  lines 
Lq{x^.  This  takes  linear  time.  For  each  r  ^  R,  determine  in  C>(log  m)  =  0(log  n) 
time  the  index  i,  1  ^  /  ^  m,  such  that  r  and  />,•  are  in  the  same  strip.  Then  we  can  in 
0(1)  time  determine  the  critical  point  y  i  Lq  corresponding  to  the  {Lq,  Pi,  .  .  .  ,Pm)- 
path  that  passes  through  r.  Also  for  each  z  =  1,  .  .  .  ,m  — 1,  we  can  determine  the 
critical  point  corresponding  to  the  {Lq,  pi,  .  .  .  ,p;^-path  where  (x,_i,  a:,-,  j:,+  i)  are 
colinear  in  0(1)  time,    a 

For  two  points  j:,  y  on  the  line  Lq{xi)  (for  any  z),  it  is  convenient  to  say  that  they 
are  {Lq,  pi,  .  .  .  ,p ^-equivalent  if  the  path  through  x  and  y  both  belong  to  the  same 
critical  interval.  We  define  the  critical  interval  containing  a  chain  to  be  the  one  which 
contains  the  first  node  Xq. 

We  now  describe  the  decomposition  process:  Let  C  =  (Cj,  .  .  .  ,C„)  be  the 
balanced  (Lq,  Li)-chain  below  a  set  R  of  points  with  pivots  pi,  .  .  .  ,p„  and  nodes 
{xq,  Xi,  .  .  .  ,x„),  m  >  0.  C  determines  a  corresponding  {Lq,Pi.  .  .  .  ,p„)-path, 
{xq,  Xi,  .  .  .  ,x„).  If  m  is  even  (resp.  odd),  we  consider  moving  the  point  Xg  along  Lq 
in  the  direction  towards  (resp.  away  from)  the  intersection  of  Lq  and  Lj.  We  move  Xq 
until  it  reach  the  first  critical  point  Xq. 
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Figure  13.  [xi^,  Xjt+j]  passes  through  r 


Recall  that  if  (x'q,  x[,  .  .  .  ,x„)  is  the  corresponding  path,  then  for  some  k  (see  Figure 
13) 

(a)  either  [x^,  Xi^+i]  passes  through  a  point  r  in  R  —  {pi,  .  .  .  ,p„}, 

(b)  or  Xjt,  Xjfc+i,  jrt+2  ^''^   colinear. 

It  should  be  noted  at  this  moment  that  our  decision  to  move  Xq  in  the  chosen 
direction  (depending  on  the  parity  of  m)  implies  that  x„  lies  outside  the  quadrant 
Lq  DLi    and  more  importantly,  the  segment  [x„-i,  x„]  does  not  become  parallel  to  Lj. 

First  consider  possibility  (a).  For  clarity,  we  will  initially  assume  that  possibility  (b) 
does  not  occur  simultaneously  and  that  r  is  unique.  Let  L2  be  the  directed  line  from 
Pk+l  to  r.  For  /  =  I,  .  .  .  ,k,  let  C,-  denote  the  chord  determined  by  [x,_i,  x,]. 
Observe  that 

C  =  (Cj cj 

is  a  balanced  (Lq,  L2)-chain.  However  (C^+j,  .  .  .  .C„)  does  not  represent  a  balanced 
(Lj,  Z-i)-chain  because  the  pivot  p„  is  at  the  midpoint  of  [x„-i,  x„]  but  x„  is  not  on 
Lj.  To  obtain  a  balanced  (L2,  Li)-chain,  we  continue  as  follows:  For 
/  =  ^  +  1 ,  .  .  .  ,m,  define  the  point 

Xi  =  Li(Xi)nL2(Xi) 
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Ll  =  ^2(^t) 


^2(^Jt+l) 


^1    =   ^l(-^m) 


Figure  14. 


Note  that  the  set  of  triangles 

L  XiXi  Xi  (i  =  k+l,  .  .  .  ,m) 

are  congruent.    Let  C,  be  the  chord  determined  by  [j:,_i,  x,].    Then  we  note  that 

C"    -   (Cjt+1,   .   .   .   ,Cm) 

is  a  balanced  (L2,  Li)-chain.  Consider  the  (L2,  Pk+i>  Pk  +  2'  ■  ■  ■  •Pm)-P^^^ 
corresponding  to  C" :  It  is  important  to  see  from  our  construction  that  JCj^+i  and  JC^^+i  ^'^^ 
in  the  same  critical  interval  with  respect  to  (L2,  Pt+i.  Pk  +  2'  ■  ■  ■  -Pm)'  '•^-  ^^  ^^  move 
from  Xi^+1  to  Xjt  +  i,  the  (L2,  Pt  +  i,  Pk  +  2'  ■  ■  ■  'Pm)'P^^^  encountered  along  the  way  are 
non-critical.    And  also  C   remain  unchanged  and  balanced. 

Now  consider  possibility  (b)  where  again  we  initially  assume  for  simplicity  that  the 
k  such  that  Xj^,  JCj^+j,  Xi^  +  2  are  colinear  is  unique.  Let  L2  be  the  line  through  x^  and 
X/^+l.  As  before  we  immediately  obtain  an  (Lq,  L2)-chain  C  of  length  k.  It  is  not  hard 
to  see  that  we  can  define  an  (L2,  Li)-chain  C"  of  length  m  —  k  —  2  by  the  same  method  as 
above. 


23 


This  completes  our  decomposition  process  for  C.  The  resulting  pair  of  chains  C 
and  C"  will  be  called  the  decomposition  of  C-  It  should  be  noted  however,  that  C  and 
C"  are  in  general  not  below  the  set  R.  Rather,  C  and  C"  are  below  some  sets 
/?(Lo,  L2)  a°d  ^(^2.  ^1).  respectively,  where  R{Lo.  L2)U/?(L2,  Li)  =  ^.  As  it  turns 
out  in  our  application,  we  do  know  the  sets  R(Lq,  L2)  and  R(L2,  L{).  To  show  the 
dependence  of  (Lq,  Li)-chains  on  the  set  R,  we  may  also  call  it  an  (Lq,  L\,  /?)-chain. 
Thus  we  have  shown  constructively  that  every  chain  of  length  m  >  0  can  be  decomposed 
into  two  chains  of  length    <  m. 

We  have  assumed  that  k  is  unique  in  (a)  and  (b)  above.  It  is  not  hard  to  provide 
the  modification  necessary  for  the  general  case.  C  is  decomposed  into  more  than  two 
chains  if  the  k  in  (a)  and  (b)  is  not  unique. 

We  next  consider  how  the  above  process  may  be  reversed,  i.e.  given  two  balanced 
chains,  check  if  they  form  the  decomposition  of  some  chain  C  and  if  so,  construct  C. 

Let  C  be  the  balanced  (Lq,  Z,2)-chain  below  R'  and  C"  the  balanced  (L2,  Li)-chain 
below  R" .  If  C'  and  C"  form  the  decomposition  of  some  C  below  the  set  ^  =  R'UR", 
then  this  could  come  about  by  the  dccocs position  process  in  one  of  the  two  ways 
corresponding  to  possibilities  (a)  and  (b)  above. 

(A)  We  first  verify  whether  case  (a)  holds:  With  the  usual  notations  for  C  and  C"  (viz. 
C  is  an  (Lq,  L2,  /?')-chain  of  length  m' ,  C"  is  an  (L2,  Lj,  /?")-chain  of  length  m" , 
Xj  and  Xj  are  the  nodes  of  C  and  C"  respectively,  etc.),  let  L2  be  a  line  through  p 
and  q  in  R.  We  first  verify  in  constant  time  that  x„-,  p ,  q,  Xq  occur  in  that  order 
in  L2  and  that  the  two  chains  do  not  intersect.    Then  we  proceed  as  follows: 


X K- 


Figure  15. 


Let  m  =  m'  +m"  +  1  and  set 
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Pi  = 


P 

q 

Pi-m'-l 


if  J  =  1 m' 

if  z  =  m'  + 1  and  m"  =  even 
if  /  =  m'  + 1  and  m"  =  odd 

if  i  =  m'  +2,  .  .  .  ,m. 


Set 


^m'  +  l   ~  Xm'\Pm'  +  l 

Xm-  +  1  =  Z,i(j:o)nLo(^m'  +  l)- 

Let  n'  be  the  (Lq,  Pi.  •  .  .  -Pm'  +  l)"P^'^  corresponding  to  moving  from  the  node 
x'„-  +  i  to  x„'  +  i  and  II"  the  (Li,  p],  pi.  .  .  .  ,p„-)-psith  corresponding  to  moving 


Figure  16. 


Verify  thatx„'  +  i  andj:„.  +  i  are(Lo,  Pi p„'  +  i)-equivalent.    and  j:„.  +  i  and 

Xq  are  (L^,  p^,  P2.  ■  ■  ■  ,/5^")-equivalent.  It  is  easy  to  see  that  these  two 
equivalence  conditions  hold  if  and  only  if  C  and  C"  form  a  decomposition  of  the 
(Lq,  Li)-balanced  chain  C. 

(B)     To  verify  if  case  (b)  holds,  we  proceed  in  substantially  the  same  way:  first  verify 
that  m"  is  even.   Let  m  =  m' +m" +  2  and 


Pi 
P 

Pi-m'-Z 


if  /  =  1 m' 

\fi  =  m'  +  l 
if  J  =  m'+2 
if  i>  m'+2 
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Set  x„-  =  x„'\p„-  +  i  and  x„' +2  =  x„' +2\p'm- +2-  We  should  verify  that  the  C  and  C" 
do  not  intersect  and 

If  I  H  t 

^m''  Pm'+l'  Pm'+2'  ^Q'  ^m'+2 
occur  in  that  order  on  L2.   This  being  the  case,  set 

We  must  verify  that  x„-+2  and  x„'+2  "^  (^0.  Pl>  •  ■  ■  ./'m') -equivalent  and  x„- +2  ^nd 
Xq  are  (Li,  p„',  p„-^i,  .  .  .  .p^) -equivalent  to  ensure  the  validity  of  the  composition 
process. 

If  C  and  C"  form  the  decomposition  of  C  then  we  can  easily  reconstruct  C.    Hence: 

Lemma  8.  Given  the  balanced  (Lq,  L2,  i?')-chain  C ,  the  balanced  (L2,  Lj,  ^")- 
chain  C"  and  the  critical  intervals  of  these  chaias,  it  takes  0(1)  time  to  verify  if  the 
balanced  {Lq,  Lj,  ^' U/?")-chain  C  decomposes  into  C  and  C" .  It  takes  0{n)  time  to 
construct  C  when  it  exists. 

Note  that  when  computing  critical  iatei-vais  for  (Lq,  p\,  .  .  .  ,/?;„) -path  on  Lq,  we 
need  to  extend  the  path  to  include  one  of  the  pivots  en  L\.  The  reason  is  if  the  chain  is 
to  compose  with  another  chain  to  its  right,  then  those  critical  intervals  are  needed. 

6.    A  Potato-Peeling  Algorithm 

By  a  potato  of  P,  we  mean  any  maximum  area  convex  polygon  contained  in  P . 
Although  our  goal  is  to  compute  the  potato  itself,  it  is  convenient  to  describe  the 
algorithm  for  computing  the  area  of  the  potato.  It  is  easy  to  modify  the  algorithm  to 
compute  the  potato  in  addition  to  its  area. 

We  now  define  some  useful  notations  and  data  structures.  For  points  r,  s  on  the 
boundary  of  P ,  let  P[r,  5]  be  the  polygon  whose  corners  arc  r  and  s  together  with  the 
corners  in  P  occuring  clockwise  from  r  to  s.  Note  that  P[r,  s]  is  a  simple  polygon  iff 
the  segment  [r,  s]  is  contained  in  P .  Let  C,-  ,•  denote  the  oriented  chord  from  v,-  to  v,  if 
the  segment  [v,-,  v^]  C  P,  otherwise  C,-  ,  is  undefined.  For  each  corner  v,-  of  P ,  define 
H,  to  be  the  set  of  extremal  chords  from  v,-: 

H,-  =  {  Cij  \Cij  is  defined  } 

Clearly  |H,|  <  n.  It  is  easy  to  compute  all  the  sets  H,  in  0(n  log  n)  time.  Let  H 
denote  the  union  of  all  the  H,'s. 

Let  {Cij,  C jc  j)  be  a  pair  of  directed  chords  from  H.  Then  define  R{C ^j,  C^_  ;) 
to  be  the  corners  of  P  clockwise  from  v,  to  v^  (inclusive).  Let  Lq  and  Lj  be  the  directed 
lines  obtained  by  extending  C ij  and  C ^.  i  respectively.    Then  we  define  a  {C ^j,  C^.   ;)- 
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chain  to  be  an  (Lq,  L-^,  R(Cij,  C^^  j))-cham  below  the  set  R(Cij,  Cjc^  /). 

Definition.  A  pair  of  chords  (C ij,  C^^  /)  is  admissible  if  the  corners  v,-,  v^-,  Vj^,  V; 
occur  in  this  cyclic  order  on  the  boundary  of  P.  A  (Cij,  Cjt  /)-chain  is  also  called 
admissible  if  (Cj  ,•,  Cj^  /)  is  admissible  and  the  first  and  last  nodes,  Xq  and  x„  (for  some 
m),  lie  in  C :j  and  C^j  respectively.  A  (C,  C")-chain  C  and  a  (C",  C')-chain  C  are 
compatible  if 

(1)  Both  C  and  C  are  admissible. 

(2)  The  last  node  x„  (for  some  m)  of  C,  the  first  node  xq  of  C  and  the  pivots  (p ,  q) 
of  C"  occur  in  the  following  order: 

Xm>  P'  ?.  ^0 

For  any  chord  Cj  ;  €  H,  we  call  the  point  mid-way  between  the  double-pivot  v,-  and  v, 
the  reference  point  of  C i  j. 

Before  giving  an  jdgorithm  for  the  potato-peeling  problem,  we  first  introduce  an 
area  measure  for  polygons.  Let  P  =  (vq,  v\,  .  .  .  ,v„_i)  be  a  polygon  that  is  not 
necessarily  simple,  with  n  corners  where  (a:,,  yj)  is  the  coordinate  of  corner 
V,-,  I  =  0,  ....  n  —  1 .    The  area  of  P  is  defined  as  follows: 

l^AREAI^P)  =   X  ^/(y.  +  l-^i) 
1  =  0 

This  is  called  the  signed  area  [Lop63]  of  P.  If  P  is  simple  polygon,  then  this  definition 
gives  the  expected  notion  of  area,  with  a  positive  sign  if  the  corners  are  given  in 
counterclockwise  order  and  negative  otherwise.  Recall  that  we  previously  define  the 
'area'  of  a  chain  or  path.  We  now  redefine  the  area  of  a  (C,  C')-chain  to  be 
AREA{{c ,  Xq,  Xi,  .  .  .  ,x„,  c'))  where  c  and  c'  are  the  respective  reference  points  of 
C  and  C'  and  the  j:,'s  are  nodes  of  the  (C,  C')-chain. 

The  following  algorithm  is  described  in  three  main  steps. 

Step  I.  We  introduce  the  matrix  A  indexed  by  pairs  of  chords  such  that  for 
C,  C  €  H,  A{C ,  C)  is  the  area  of  the  unique  balanced  (C,  C')-chain,  if  it  is 
admissible.  Otherwise  A(C,  C)  =  -».  Initially  set  A(C,  C)  to  -oo  for  all  (C,  C). 
Then  for  each  admissible  pair  (C,  C),  find  R{C ,  C).  It  is  easy  by  a  brute-force 
method  to  compute  all  the  (Cq,  Ci)-chains  of  lengths  0  or  1  in  time  0{n  ).  Note  that 
the  chain  has  length  0  precisely  where  R{C ,  OnC^nC'"^  is  empty. 

To  compute  admissible  chains  of  all  lengths,  we  proceed  in  stages.  The  first  stage  is 
the  computation  of  admissible  chains  of  lengths  0  and  1.  At  stage  i+  1  we  compute  more 
admissible  chains  by  composing  admissible  chains  computed  in  the  previous  stages.  To 
compute    a    (C,  C) -chain    at    stage    z'+l,    we    iterate    through    all    chords    C"    checking 
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whether  the  (C,  C")-  and  (C",  C')-chains  have  been  computed  in  previous  stages  and 
whether  they  form  a  decomposition  of  the  (C,  C')-chain.  If  so,  construct  the  (C,  C')- 
chain.  At  the  same  time,  we  should  compute  the  critical  intervals  for  the  (C,  C')-chain 
with  respect  to  the  enlarged  set  of  points,  RiC,  C):  this  takes  0(n\ogn)  time  as  shown 
in  Lemma  7.  In  order  to  facilitate  composition  of  chains,  we  can  determine  at  the  same 
time  the  critical  interval  to  which  the  first  node  of  the  (C,  C')-chain  belongs. 

To  analyze  the  complexity  of  this  procedure,  we  divide  the  cost  into  two  parts:  (i) 
the  cost  for  verifying  the  possibility  of  composition  and  (ii)  the  cost  of  actually  composing 
the  chains  and  computing  the  critical  intervals.  There  are  0(n)  stages.  At  each  stage, 
we  go  through  all  triples  (C,  C" ,  C)  verifying  if  it  is  possible  to  compose  the  (C,  C")- 
and  (C",  C')-chains.  Hence,  at  each  stage,  there  are  (n  j  instances  of  testing  for  a 
possible  composition  where  each  test  takes  0(1)  time,  Since  there  are  n  stages,  the  cost 
of  part  (i)  is  0(n  ).  Because  of  the  uniqueness  of  the  balanced  chains,  there  are  at  most 
0{n  )  chains  that  are  composed  during  the  entire  procedure.  The  cost  of  part  (ii)  is 
0{n  logrt),  since  each  composition  and  computation  of  critical  intervals  take  0(/zlogn) 
time.   Thus  the  total  cost  of  computing  the  A  matrix  is  0(n  ). 

Step  II.  First  we  introduce  the  matrix  M  with  entries  indexed  by  pairs  (C,  C)  of 
extremal  chords.  M{C ,  C)  is  the  maximum  area  of  admissible  (C,  C')-chains; 
M{C ,  C)  =  — 00  if  (C,  C)  is  not  an  admissible  pair  of  chords.  M  can  be  computed  in 
stages  where  in  stage  5  (^  =   1,  2,  ...,«  — 1)  we  compute  the  entries  M(C ,  C)  where 

C  ^  E,-,  C  €  Sj  and  j  —  i  =  s  (mod  n) 

Note  that  in  stage  j  =  1,  we  have  j  =  /+!  (mod  n)  and  the  constraint  that  (C,  C) 
forms  an  admissible  pair  implies  that  C  must  be  the  chord  C,  ..  Let  c  and  c'  be  the 
reference  points  of  C  and  C  respectively.  Then  M(C ,  C)  is  given  by  the  area  of 
A  cvjc' . 


Figure  17. 
In  general,  for  stage  s  >  1,  we  use  the  recursive  formula: 
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M(C,C')  =  msixi  A(C .  C),  maxi  M(C .  C")  +  M(C" ,  C)  +AREA(^cc"c')  }  }, 

where  C"  ranges  over  "S./^,  k  has  the  range  i  <  k  <  j  and  c"  is  the  reference  point  of 
C" .  (See  Figure  18.)  To  justify  this  formula,  note  that  A/(C,  C)  is  either  determined 
by  the  balanced  (C,  C) -chain  or  else  it  is  determined  by  a  chain  that  has  an  extremal 
chord  C"  ^  Hj. 


A  cc  c 


Figure  18. 

Note  that  we  do  not  check  compatibility  between  M(C ,  C")  and  M(C" ,  C).  If  they  are 
not  compatible,  then  they  cannot  form  an  optimal  chain.  Thus  it  is  not  necessary  to 
exclude  them  from  the  maximization.  It  takes  0{n  )  steps  to  carry  out  the  maximization 
for  each  entry  of  M .   Thus  the  whole  matrix  M  takes  0{n  )  steps  in  total  to  compute. 

Step  in.  The  potato-peeling  problem  is  now  solved:  We  know  that  there  are 
m  '^  2  extremal  chords  that  form  part  of  the  boundary  of  the  potato.  It  can  be  found  in 
0{n  )  time  as  follows: 

max{A/(C,  C')+M{C' ,  C)  :  C  ^  H,-,  C  €  Hy,  for  all  /and;}. 

Notice  again  we  do  not  check  compatibility  between  M{C ,  C)  and  M(C' ,  C). 

7.    Potato  Peeling  -  Perimeter  Measure 

The  potato-peeling  problem  under  the  perimeter  measure  can  be  solved  in  the 
slightly  better  time  bound  0(n  )  by  using  essentially  the  same  techniques.  However, 
some  additional  properties  here  make  a  much  simpler  algorithm  possible. 

First  of  all,  we  give  a  finiteness  criterion  for  the  problem.  Let  us  consider  the 
simplest  case  of  the  perimeter  optimization  problem  where  the  given  polygon  has  just  one 
reflex  corner,  say  vq.  As  with  the  area  measure,  the  problem  is  to  determine  the  chord 
C  =  [c,  c']  through  Vq  that  maximizes  the  perimeter  of  the  convex  polygon  PClC 

We  shall  prove  that  C  may  be  assumed  to  be  one  of  the  two  extremal  chords  of  the 
butterfly  containing  C.    Thus  in  contrast  to  the  area  measure,  we  need  not  consider  a 
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third  possibility  (such  as  C  being  balanced). 


'  ^^^ra~ 


Figure  19.  The  perimeter  determined  by  a  chord  in  a  butterfly 

Notation:  Given  any  two  points  a  and  b  on  the  boundary  of  P ,  let  S[a,  b]  denote 
the  length  of  polygonal  path  clockwise  between  a  and  b  if  the  path  contains  no  reflex 
vertex,  and  —  <»  otherwise. 

Consider  the  butterfly  B  that  contains  C .  Let  Lq  and  L\  be  the  two  supporting  lines 
of  S.  If  Lq  and  Lj,  the  two  supporting  line  of  B  arc  parallel  then  the  perimeter  is 
independent  of  the  choice  of  C;  we  are  done.  Hence  assume  that  Lq  and  Lj  intersect  at  a 
point  o.    Referring  to  Figure  19,  let 

a  =  Z.  a'ovQ, 
P  =  TT-Z.  aOVQ, 


9  =  Z.  ovqc. 

Note  that  0  ^  a  <  ^  <  it.    Then  9  satisfies  a  <  9  <  (3.    And  by  the  law  of  sines  we 
have 


\oc\ 


|voc| 


\ovq\ 


sin  9         sin  (ir-p)         sin  (p-9)  ' 

\vqc'\    _  |oc'| IovqI 

sin  ct         sin  (tt  — 9)         sin  (9-a) 
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Without  loss  of  generality,  let  us  assume  |ovo|  =  1.   Hence  we  have 

I  M       I       Ij.1      'I  sin  p  sin  g 

\cc  I  =  |cvo|+  Ivoc  I  =    3in(3l9)  +  sin(9-a)' 

II  sin  9 
\oc\  = 


\oc-    = 


sin  (3-9)  ' 

sin  9 
sin  (9-a)  ' 


We  now  obtain  the  perimeter  L  as  a  function  of  9.   It  suffices  to  show  that  L  has  no 
local  maximum  to  prove  that  C  is  extremal.    We  consider  the  following  two  cases: 

(1)  B  is  an  V-buttcrfly: 

C*"   contains  o  and  the  perimeter  L(9)  of  the  convex  polygon  determined  by  the 
chord  C  is 

L(9)  =  \ac\'^\cc'\+\c'b'\  +  S[b' ,a] 

=  \oc\+  \cc'\+\c'o\-  \oa\-  \b'o\  +  S[b',a] 

=  |oc|+ |cc'|+ |c'o|  +  constant 

sin  9  +  sin  a   ,   sin  9  +  sin  B   , 

=  — : — —  +  — — r^ — r-^  +  constant . 

sin  t,9-a)  sin  (p-9) 

dL  _   cos  9  sin  (9  — a)  — (sin  9  +  sin  a)  cos  (9  — a) 
dQ  sin2(9-a) 

cos  9  sin  (p-9)  +  (sin  9  4- sin  p)  cos  (p-9) 
sin^(3-9) 

sin  g  (14- cos  (9-a))       sin  P  (1 +cos  (p-9)) 
sin2(9-a)  sin2(p-9) 

_  -sin  g sin  P 

l-cos(e-a)        l-cos(p-9)' 

d^L    _    sin  g  sin  (9  — g)         sin  P  sin  (p  — 9) 
dQ^  (1-cos  (9-g))2       (1 -cos  (p-9))2  " 

d^L 

Since  g,  P,  9  — g,  and  p— 9  arc  all  in  the  first  two  quadrants,  r"  —  0- 

dQ 

(2)  B  is  a  A-butterfly: 
Similiarly  we  have 
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L(9)    =  \a'c'\+\c'c\+\cb\  +  S[b,a'] 

=  -  \c'o\+  \c'c\-  \co\+\a'o\+  \bo\  +  S[b,  a'] 

—  —  |oc' I  +  |c'c|— |co  I  +  constant 

sin  a  — sin  9   ,   sin  B  — sin  8   , 

=  — : — — : — I .    ,^ — — — h  constant, 

sin  (e-a)  sm(3-e) 

dL  _  sin  g sin  p 

dQ   ~    l  +  cos(9-a)       l  +  cos(3-8)' 

d^L  __    sin  a  sin  (9- a)        sin  P  sin  (^-9)    ^ 
dQ^  ~   (l  +  cos  (9~a))2       (1 +  cos  ((5-9))^ 

In  both  cases  Z-(9)  does  not  have  any  local  maxima.  We  have  proved  the  following 
lemma. 

Lemma  9.  For  a  polygon  with  one  reflex  corner,  the  maximal  perimeter  is 
determined  by  an  extremal  chord  through  the  reflex  corner. 

A  chord  is  semi-extremal  if  it  shares  a  common  endpoint  with  an  extremal  chord. 

Lemma  10.  Let  C  =  (Cj,  C2,  ■  ■  •  .C„)  be  a  optimal  chain  of  chords  for  a  series 
of  butterflies  (Sj,  B2,  .  .  ■  .B„).  For  each  /,  1  ^  j  <  m,  either  both  C,-  and  C,  +  i  are 
extremal,  or  one  of  them  is  extremal  and  the  other  is  semi-extremal. 

Proof.  By  definition  of  a  chain,  C,  and  C,+  i  must  intersect.  We  say  two  line 
segments  overlap  if  their  intersection  has  positive  length.    There  are  two  cases: 

(1)     B  j  and  S,  +  i  'intersect  fully'  i.e.  the  forward  tip  of  5,  does  not  overlap  the  backward 
tip  of  5,+ 1. 
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Figure  20.  TWQ  butterflies  that  intersect  fully 

If  we  consider  C,-  to  be  fixed,  then  Lemma  9  implies  that  C,+i  must  be  extremal  in 
PDC^ .  But  since  Bi  and  S,  +  i  do  not  share  a  tip,  we  see  that  C,+  i  must  in  fact  be 
extremal  in  P.   Similarly,  C,-  is  extremal. 

(2)     The  forward  tip  of  5,-  overlaps  the  backward  tip  of  S,  +  i. 

First  we  assume  that  the  two  overlapping  tips  are  identical  (i.e.  [a,  b]  in  Figure  21). 
Let  c  be  the  intersection  of  C,-  and  C,  +  i.  It  is  not  hard  to  see  that  Lemma  9  implies 
that  c  lies  on  the  boundary  of  B,ni5,+  i.  Suppose  c  does  not  lie  on  the  shared  tip  T 
{T  =  [a,  b]).  Again  Lemma  9  implies  that  C,-  and  C,  +  i  must  be  extremal  in  their 
respective  butterflies. 


Figure  21.  Two  butterflies  with  overlapping  tips 


Now   consider  the  case  where  c   is  on  the  shared  tip.    The  perimeter  L  =  L(c) 
determined  by  c  is 
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=  (|cc,|+|cia,|+|cZ;|)  +  (|cc,+  i|+h+iZ7i+i|+|ca|) 

=  L,<c)+L,+i(c)  + constant, 
where  L,<c)  and  L,  +  i(c)   are  perimeter  functions  for  the  butterflies  S,-  and  fi,  +  i 

respectively.    From  the  proof  of  Lemma  9  we  know  that -r-  ^  0  and r-  S  0. 

So  L(c)  has  positive  second  derivative  and  is  unimodal.  We  conclude  that  L  has  no 
local  maximum  for  c  in  the  range  [a,  b]  and  the  maximal  perimeter  is  determined 
by  an  endpoint  of  the  shared  tip.  Hence  both  C,-  and  C,  +  i  are  extremal.  Finally  if 
we  drop  the  assumption  that  the  overlapping  tips  of  fl,-  and  S,  +  i  are  identical,  the 
analysis  can  be  modified  in  the  obvious  way  to  show  that  one  of  C,  and  C,  +  i  is 
extremal  smd  the  other  is  semi-extremal. 

We  conclude  that  in  all  cases  at  least  one  of  C,-  and  C,  +  i  is  extremal  and  the  other 
is  either  extremal  or  semi-extremal.    □ 

Notation:  For  vertices  v,-,  v  of  P,  let  P[i,  j]  denote  the  polygon  formed  from  the 
vertices  of  P  clockwisely  between  v,-  and  Vj.  P[i,  ;]  is  undefined  if  [v,-,  vj\  is  not  fully 
contained  in  P.  For  chords  C,-  ^  H,-,  Cj  i  Hj,  let  P[i,j,  C i,  C j\  denote  the  connected 
polygonal  region  in  P[/,  yjOC,^  HC/"  that  is  bounded  by  [v,-,  v.J.  Note  that 
P[i,  i^f^C^  r\Cj'  need  not  be  a  connected  region. 

Define  A/°j(C,-,  C ,)  to  be  the  length  of  the  polygon  P[i,j,  C,-,  Cj]  (not  counting 
the  edge  [v-,  vJ)  if  P[i,  j]  is  defined;  and  — oo  otherwise.  M  can  be  computed  in 
0(n  )  time.  Let  M'jiCi,  Cj)  denote  the  perimeter  of  the  largest  convex  polygon 
contained  in  P[i,j,  C,-,  C,]  (again  the  edge  [j ,  i]  is  not  counted).  For  fixed  /'  and  y,  we 
can  regard  each  M ij(C i,  Cj)  as  an  n  X  n  matrix. 

Now  we  arc  ready  to  present  an  0(n  )  time  algorithm  for  finding  the  largest 
perimeter  potato.  Again,  instead  of  computing  the  potato  itself,  we  will  pretend  that  we 
are  only  computing  its  perimeter. 

Our  problem  is  essentially  reduced  to  computing  M* .  To  show  this  reduction, 
suppose  that  the  potato  in  P  is  determined  by  a  set  of  chords  where  at  least  two  are 
extremal.    So  the  maximum  perimeter  is  given  by 

max  {  m: ,(C,..  Cj)  +  Mj  ,(C,.  C,)  }. 

where  U  is  the  set  of  reflex  vertices  of  P.    This  expression  can  be  evaluated  in  0(n  ) 
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time,  given  M*.  The  case  where  there  is  at  most  one  extremal  chord  in  the  potato  can  be 
done  in  0{n  )  time  using  a  brute  force  method.  Note  that  this  case  implies  that  the 
potato  is  determined  by  at  most  three  chords.   It  remains  to  show  how  to  determine  A/*. 

To  compute  A/*,  we  define  two  additional  nXn  matrixes  M  and  A/  .  With 
/,  j,  Ci,  Cj  as  before,  and  for  m  =  1,2,  we  define  M^jiC,,  Cj)  to  be  the  perimeter  of 
the  largest  convex  polygon  contained  in  P[i,j,  C,-,  Cj]  determined  by  at  most  m  semi- 
extremal  chords  in  addition  to  C,,  C,  (but  no  other  extremal  chords).  Again  the  length 
of  the  edge  [i,  j]  is  not  counted  as  part  of  the  perimeter.  We  then  have  the  following 
formula: 

Mlj{Ci,Cj)=         max    ^    {  A/J/C,-,  C,),   A/.^/C,,  C,),   A/;.  i(C,-.  C^) +M;,/Ci,  C,)  }. 

This  formula  is  justified  by  the  fact  that  if  the  potato  in  P[i,  j,  C,-,  C,]  is  determined  by  3 
or  more  chords  then  at  least  one  is  extremal.  As  in  the  case  of  the  area  measure,  we  can 
recursively  compute  the  entries  of  A/*  in  n  stages.  We  initialize  A/*  to  M  ,  i.e. 
M*j(Ci,  Cj)  =  M ij{C i,  Cj),  to  start  this  recursion.  It  is  easy  to  see  that  0(n  j  time 
suffices  for  the  overall  computation,  assuming  the  availability  of  M  ,  M    and  M  . 


Figure  22.  Computing  M 

Both  M  and  M  can  be  computed  in  0{n  )  steps.  We  now  describe  briefly  how 
M  is  computed.  For  reflex  vertices  v,,  Vj  and  extremal  chord  C,  (:  s,,  C,  d  S,,  let  C 
and  C'  be  the  two  extremal  chord  that  determine  MfjiCi,  Cj).  We  observe  that  C  and 
C  are  disjoint  and  share  endpoints  respectively  with  C,  and  C,.  (See  Figure  22.)  It  is 
clear  that  C  and  C  can  be  found  in  0(n  )  steps  by  an  exhausted  search,  provided  that 
all  the  semi-extremal  chords  are  precomputed.    Since  there  are  0(n  )  entries  in  M  ,  it 

ft  0  1 

take  0(n  )  steps  to  compute  M  .    We  can  similarly  compute  A/  . 


8.    Conclusion 

This  paper  gives  the  first  polynomial-time  solution  to  the  potato-peeling  problem. 
We  have  introduced  the  interesting  geometric  concept  of  balanced  chains  which  holds  the 
key  to  the  problem.    Computing  these  chains  is  the  bottleneck  to  a  faster  algorithm  for 
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the  problem.    Our  solution  also  exploits  dynamic  programming  in  several  key  steps. 
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